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ABSTRACT 

Gravitational lensing can provide pure geometric tests of the structure of space-time, for instance by 
determining empirically the angular diameter distance-redshift relation. This geometric test has been 
demonstrated several times using massive clusters which produce a large lensing signal. In this case, 
matter at a single redshift dominates the lensing signal, so the analysis is straightforward. It is less 
clear how weaker signals from multiple sources at different redshifts can be stacked to demonstrate 
the geometric dependence. We introduce a simple measure of relative shear which for flat cosmologies 
separates the effect of lens and source positions into multiplicative terms, allowing signals from many 
different source-lens pairs to be combined. Applying this technique to a sample of groups and low-mass 
clusters in the COSMOS survey, we detect a clear variation of shear with distance behind the lens. This 
represents the first detection of the geometric effect using weak lensing by multiple, low-mass systems. 
The variation of distance with redshift is measured with sufficient precision to constrain the equation 
of state of the universe under the assumption of flatness, equivalent to a detection of a dark energy 
component fix at greater than 99% confidence for an equation-of-state parameter —2.5 < w < —0.1. 
For the case w — —1, we find a value for the cosmological constant density parameter JIa — O.SS^Q ^g'' 
(68% C.L.), and detect cosmic acceleration (go < 0) at the 98% C.L.. We consider the systematic 
uncertainties associated with this technique and discuss the prospects for applying it in forthcoming 
weak- lensing surveys. 

Subject headings: cosmology: observations - gravitational lensing - dark matter - large-scale structure 
of universe - galaxies: groups: general 
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1. INTRODUCTION 
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The current evidence for a do minant dark e nergy 
component in the universe (e.g. iPercival et al.l 120101 : 
iLarson et al.ll20lit ISullivan et al.ll201lD leaves cosmoloev 
in a uncomfortable situation. Given the concurrent evi- 
dence for cold dark matter (CDM) and an additional field 
driving inflation, it seems several radical new components 
of physics are required to explain the present-day state 
of the universe, with little detailed observational infor- 
mation so far as to their precise nature. Clarifying the 
nature of dark energy is particularly challenging. Ob- 
servationally, the effect of dark energy on the equation 
of state is very close to that of a cosmological constant 
A. To determine anything else about this component 
requires very precise tests and a rigorous elimination of 
systematics. 

The simplest evidence for dark energy comes from 
measurements of the geometry of space-time, or equiv- 
alently dist ance as a function o f redshift, either from 
supernovae (|Sullivan et al.l |20 111) , which measure lumi- 
nosity distance , or baryon acoustic oscill ations (BAO) 
(jPercival et al.l 120101: iBeutler et all 120111) or _the_cos- 
mic microwave background (CMB) (|L arson "etaI1[20Tll) . 
which measure angular diameter distance. Gravita- 
tional lensing also provides tests of cosmology, measuring 
both the matter distribution and how gravity deflects 
light on large scales. There has been much emphasis 
on c osmological lens i ng tests using c osmic shear (e.g. 
■Massey et atl l2007at iFue^l.l 120081: iSchrabback et al.l 
120101 ISembdoiir et al.l 120111 ). which is sensitive both 
to the matter distribution and to space-time geometry. 
While the theory of CDM structure formation makes 
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fairly robust predictions as to the matter distribution, 
and thus the lensing potential, on scales larger than 
galaxies, it also assumes additional physics associated 
with inflation, such as Gaussian distribution of initial 
fluctuations with an almost-scale invariant power spec- 
trum. By separating out the geometric effects of lens- 
ing from the properties of the density field, one can in 
principle derive more general constraints on geometry, 
independent of these assumptions. 

The underlying idea of geometric lensing tests is 
straightforward. First, the strength of lensing is mea- 
sured behind an object as a function of redshift. Then, 
taking ratios of the lensing signal at different redshifts, 
the dependence on redshift gives a measure of space-time 
geometry via the angular diameter distance-redshift re- 
lation. The attraction of this technique is that (at least 
in principle) any uncertainties in the exact form of the 
lens potential cancel out. A cosmological constant (or 
dark energy with a similar equation of state) manifests 
itself by increasing the distance to an object at a given 
redshift. In geometric lensing tests, the amplitude of the 
lensing distortion provides an independent estimate of 
the distance that the source lies behind the lens. Thus 
by measuring average distortion versus offset in redshift, 
one can constrain the value of A or more generally fix, 
the dark energy density parameter. 

In practice, the specific implementations of this idea 
fall into two broad categories. The most developed tests 
use one or a few massive clusters to provide the lensing 
signal, and measure the position of strongly lensed arcs 
to determine the strength of lensing as a function of red- 
shift behind the cluster. Th i s met hod was first discussed 
in detail bv iLink fc Pierce! Jl998i) (although the idea is 
much older, e.g. Rcfsdal 196^, who assumed the simplest 
singular isothermal sphere (SIS) potential for the cluster. 
The m ethod was subsequently revisited by Golse et al. 
(j2002l ) , who considered several sources of systematic er- 
ror, in particular the effect of substructure and irregular 
cluster mass distributions. They applied the test to the 
galaxy clusters AC114 and Abell 2218, obtaining con- 
straints < < 0.33 on the matter density parameter 
and w < —0.85 on the dark energy equation-of-state pa- 
rameter JSoucail et alll2004[ ). The method was also used 
bv lSerenor(|2002l ) on the cluster CL 0024-^1654, providing 
evidence fo r a flat accelerating cosmology. Most recently, 
I Julio et all ()2010l) have derived tight constraints using 28 
images from 12 multiple image families in Abell 1689. 

A related method uses weak lensing to measure the 
amplitude of the lensing sig nal as a function of red- 
shift behind a clus ter (see iLombardi fc BertinI 119991 : 
Gautret et al.l l2000L for early forms of this test). 



Wittman et al. r iooi first used the weak lensing signal 



to determine the redshift to a cluste r by lensing alone, 
and it has subsequently been used by lGavazzi fc Soucaill 
(|2007( ) to estimate redshifts to a dozen clu sters in the 
CFHT LS Deep fields. Most recently Mede zinski et al.l 
()2011t ). have demonstrated the effect for 3 massive clus- 
ters using rough photometric redshifts (photo-zs) to de- 
termine mean redshifts to different samples of lensed 
galaxies. This work should provide interesting con- 
straints on dark energy when extended to their full sam- 
ple of 25 clusters. 

A second category of tests uses the combined signal 
from many less massive halos as the source of the lens- 



ing effect. iJain fc Tavloil ()2003[ ) proposed the first such 
test, the 'cross-correlation tomography' method. Here 
the mass distribution is inferred statistically from the 
foreground galaxy distribution, while the lensing signal 
is measured with weak shear in two background samples; 
ratios of the galaxy-shear cross-correlation functions for 
the two samples then probe ge ometry. Variants of this 
met hod have been devel op ed bvlBerns tein fc Jain' (2004') 
and IZhang et al.l (|2005| ). iTavlor et al . (2007) proposed 
applying this technique behind clusters using ratios of 
individual shear measurements, rather than correlation 
functions. Their revised met hod was demonstrate d to- 
gether with 3D cosmic shear in lKitching et al.l (|2007( ) and 
sy stematics and error fo recasts were discussed in detail 
in lKitching et al.l dlOOS"). 

The two sorts of tests have different advantages and 
disadvantages. Tests using a small number of mas- 
sive clusters benefit from significant lensing signal and 
external constraints on the form of the mass distribu- 
tion from optical or X-ray data, but suffer from cos- 
mic sample variance and uncertainties in the cluster pro- 
file. Since these tests probe only one or a few lines of 
sight, real structures in front of or behind the cluster - 
voids or other halos - will enhance or reduce the sig- 
nal at particular redshifts. Although ideas for model- 
ing t he line of sight are emerging (e.g. Hoekstra et al] 
1201 It ), the results are not yet satisfactory, and the cur- 
rent practical consensus is to average the signal com- 
ing from different clusters. Simulations suggest that at 
least ~10 massive clusters with many multiple-image 
systems each are needed to overcome cosmic sample 
variance (iDalal et al.1l2005t iGilmore fc Nataraianll200l 



iD'Aloisio fc NataraianI l2011i r~ Furthermore, massive 
clusters are relatively rare, so only a fraction of the to- 
tal lensing signal from all cosmic structure can be used. 
Tests using a larger number of less powerful lenses sam- 
ple more of the total lensing signal but require large 
survey areas with accurate photometric redshifts, and 
may be subject to more subtle uncertainties in the fore- 
ground mass distribution. Furthermore, since the geo- 
metric term in the lensing equation depends on source 
and lens redshifts separately, it is not immediately clear 
how to stack the results from large samples of lenses in a 
simple way. The analysis is thus less intuitive, making it 
harder to spot unanticipated systematics in the results. 

The COSMOS survey provides an interesting data set 
with which to explore geometric lensing tests. COS- 
MOS has an unusual combination of a high density of 
sources with accurate lensing shape measurements, and 
accurate photo-zs for a large fraction of these sources. 
Unfortunately the COSMOS field has no really massive 
clusters in it; the largest cluster has an e stimated mass 
of 2.5 X IO^Mq (jFinoguenov et al.ll2007D . 10 times less 
than the largest str ong-lensing cluster s, and is at a red- 
shift of z = 0.73 (|Guzzo et a,l.l l2007l ) where lensing is 
past the peak in sensitivity, given the source redshift 
distribution. The lensing signal in the COSMOS field 
comes i nstead from many low-significance, group-size d 
haloes (jFinoguenov et al.ll2007l: iLeauthaud et aLlboiOD . 
Collectively these systems could still produce a large lens- 
ing signal to test geometry, however, provided the signal 
could be stacked. 

Here we propose a simple method for stacking the sig- 
nal from multiple lenses into a single measure of geom- 
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etry, in effect the relation between comoving or angu- 
lar diameter distance and redshift. Applying this new 
'stacked shear ratio test' to X-ray selected groups in the 
COSMOS field, we obtain a clear detection of the ge- 
ometric signal and derive significant constraints on the 
dark energy density parameter fix- While the COSMOS 
field is probably too small to overcome sample variance 
limitations, the magnitude of our statistical errors illus- 
trates the future promise of this technique. 

The outline of paper is as follows: in section 2 we 
present the basic data, including the sample of lensing 
groups, and the source redshifts and shape measurements 
of the COSMOS lensing catalogue. In section 3 we intro- 
duce the stacking technique and discuss optimal weight- 
ing for this technique. In section 4, we use the stacked 
shear ratio test to derive parameter constraints on the 
density of dark energy fix and the equation-of-state pa- 
rameter w, and discuss possible systematics for this test. 
In section 5 we summarize our results and discuss fu- 
ture prospects for applying the stacked shear ratio test 
to other weak lensing surveys. Throughout the paper 
we refer our results t o the WMAP 7-year mean param- 
eter values derived in iLarson et al.l ()2011[ ). taking a flat 
cosmology with Qa = 0.73, = 0.27, Hq = 70/i7o 
kms-i Mpc"i (WMAP7 hereafter) as our baseline. 

2. DATA 

2.1. The COSMOS survey and group sample 

The COSMOS survey (iScoville et al.ll2007al) brings to- 
gether panchromatic imaging from X-ray to radio wave- 
lengths, including the largest contiguous area observed 
with the Hubble Space Telescope, and deep optical spec- 
troscopic observations. The field covers an area of 
1.64 square degrees centered at 10:00:28.6, -^02:12:21.0 
(J2000) and contains identified groups, clusters and 
larger structures sp anning a wide range in redshift 
(jScoville et al.ll2007H ). 

We consider the gravitational lensing signal behind a 
sample of gal axy groups selected orig inally via their X- 
ray emission ()Finoguenov et al.1l2007() and updated us- 
ing a co mbined mosaic of imaging from XMM-Ne wton 
(1.5 Ms. lHasinger et al.ll2007l:ICappelluti et al.ll2009D an d 
the Chandra observatories (1.8 Ms. lElvis et al.l l2009t ). 
Groups are detected from the combined X-ray mosaic 
using a wavelet filter, which can result in centering uncer- 
tainties of up to 32". The distribution of galaxies along 
the line of sight to each X-ray detection is searched for 
a red sequence overdensity to determine the group red- 
shift, wit h spectroscopic redshifts used for subsequent re- 
finement (jFinoguenov et al]|2007[ ). Group members are 
selected based on their photometric redshift and prox- 
imity to the X-ray centroid, using an algorithm tested 
exte nsively on mock cat alogs and spectroscopic subsam- 
ples ()George et al.ll20lTl ). Stellar masses of the member 
galaxies are determined from multiwavelength data (see 
iLeauthaud et al. 2011 for details). From an initial list 
of members, group centers are then redefined around the 
most massive group galaxy within the NFW scale ra- 
dius of the X-ray centroid (MMGGscaio), which optimizes 
the weak lensing signal at small radii (George et al. in 
prep.). For the majority of our groups this gives centers 
which agree with the X-ray centroid; a minority (approx- 
imately 20%) of groups show significant offsets between 
the most massive galaxy and the X-ray centroid. These 



offsets could be due to observational problems (such as 
low signal-to- noise in the X-ray or optical data), or they 
might indicate unrelaxed, low-concentration groups with 
poorly defined physical centers, such as recent mergers. 
We will consider below both the full set of groups and 
a 'restricted' set which excludes the systems with sig- 
nificant offsets. The centering algorithm will be dis- 
cussed further in a forthcoming paper (George et al. in 
prep.). The full X-ray group sample, together with 
derived properties, will be made available through the 
NA SA/IPAC Infrared Science Archive (IRSA) website^ 
(see lGeorge et al.ll2011l for details). 

We restrict the lens sample to groups at z < 1 to en- 
sure the reliability of X-ray detections and optical associ- 
ations, as well as good photometric redshifts for identify- 
ing members and centers. We further cut out of the sam- 
ple poor groups, groups with centroids affected by mask- 
ing, and possible mergers (this corresponds to taking only 
group s with FLAG_INCLUDE = 1 as defined in lGeorge et al.l 
l2011h . Our final sample consists of 129 systems (105 
in the restricted set) spanning a rest-frame 0.1-2.4 keV 
luminosity range between lO''^ and 10*^erg/s, with esti- 
mated virial masses of 0.8 X 10^3 - 2 X 10" /i^ M0, virial 
radii of 0.4-0.8 /ifg^Mpc, and projected angular sizes of 
l'-6'. Fig. [1] shows the mass, physical size and angular 
size for the groups in the sample (note units have been 
converted from the value Hq = 72 km s^^ Mpc~^ used 
in the catalog to Hq = 70 km s~^ Mpc~^). The virial 
radius is taken to be i?200c, the radius within which the 
mean density is equal to 200 times the critical density 
pdz) at the redshift of the group, and the virial mass 
is taken to be M200C, the mass enclosed within i?200c- 
These masses and radii are estimated from the X-ray 
data, usi ng X-ray scaling re lations c alibrated with lens- 
ing data (iFiiioguenov et al.l[2007l : ILeauthaud et al.ll2010l : 
iGeorge et al.ll2011h . 

2.2. Weak lensing galaxy shape measurements 

High resolution imaging of the COSMOS field was 
obtained with the Hubble Space Tel escope between 
October 2003 and Ju ne 2005 ((Scoville et all l2007cl : 
iKoekemoer et all I2007T ) . The main program consisted 
of 575 slightly overlapping pointings of the Advanced 
Camera for Surveys (ACS) Wide Field Camera (WFC) 
taken with the FSIAW (approximately /-band) filter. 
At each pointing, four slightly dithered exposures of 
507 seconds were obtained. Any cosmetic defects and 
refiection ghosts were carefully masked by hand. Us- 
ing Version 2.5.0 of the SExtractor photometry package 
(|Bertin fc A rnouts 1996), in a Hot- Cold configuration on 
the stacked images, we detected compact objects in a 
0.1 5" diameter aperture d own to FSIAWab = 26.6 at 
5(7 (ILeauthaud et al.ll2007f) . 

We measure the sh apes of galaxies in th is catalogue us- 
ing the RRG method (iRhodes et a l." 20001. largely follow- 
ing the analysis pipeline of ILeauth aud e t al., (2007 ) and 
calibrated against sii nulated ACS images generat ed with 
the s image package ()Massev et al.ll2004 l2007bf) . How- 
ever, we now include two significant improvements on 
this earlier work. 

First, we correct trailing in the ACS images due 
to charge transfer inefficiency (CTI) via a physically- 

^ http:/ /irsa.ipac.caltech.cdu/Missions/cosmos.html 
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Fig. 1. — Top panel: group mass versus redshift. Middle panel: 
estimated virial radius in h^^ Mpc (filled circles). The solid squares 
show the estimated scale radius, which is typically 1/4-1/5 of the 
virial radius. Bottom panel: estimated angular extent of the virial 
radius. In each case, only the 129 groups with well-determined 
centers and redshifts are shown, as described in section 2.1. A 
WMAP7 cosmology is assumed. 

motivated readout model (jMassev et al.ll20l"o( ) that acts 
at the pixel level, rather than a parametric scheme at 
the catalog level. This moves electrons in the raw expo- 
sures back to where they should have been read out. The 
method achieves a 97% level of correction and is robust 
to variety in galaxy morphology, local galaxy density and 
sky background level. After correction, residual shears 
are well below statistical measurement precision. 

Second, we model the Hubble Space Telescope's point 
spread function (PSF) as a function of chip position, 
telescope focus offset, and velocity aberration factor 
(VAFACTOR). The latter two parameters reflect HST's 
thermal condition during each exposure: slight ex- 
pansion and contraction changes the PSF. Following 
[Rhodes et all (|2007l ). we measure the focus offset (the 
distance between the primary and secondary mirrors) 
with a precision of 1/im by comparing the apparent 
sha pes of ^ 10 s tar s in each exposur e to TinyTim mod- 
els (|Kristl 120031 ). iJee et al.l (|2007t ) found that focus 
offset correlates with the first Principal Component of 
PSF variation and accounts for 97% of the power, while 
ISchrabback eFaH f2010) found that VAFOCUS corre- 
lates with the next Principal Component. We measure 
the shapes of all stars in the COSMOS imaging, then in- 
terpolate between them using all four mea sured parame- 
ters. T his improves residuals compared to iMass ev et al.l 
(|2007al) . and we retain this physically- motivated ap- 
proach rather than relying solely on Principal Compo- 
nent Analysis. 

We have also revised our method for determining the 
variance of the tangential shear slightly. This is now 
determined empi r ically , as described in section 3.5 of 
iLeauthaud et al.l (|2011t ). Galaxies are binned by S/N 
and magnitude, and the total variance of the shear com- 
ponents 7i and 72 is measured directly in each bin. This 



empirical derivation of the shear dispersion includes both 
the scatter due to intrinsic variations in galaxy shape, 
and the additional scatter due to shape measurement 
errors. We find that the shear dispersion varies from 
cr^j, ~ 0.25 for bright galaxies with high S/N to a.y ^ 0.4 
for faint galaxies with low S/N. These measured values 
may be very slightly overestimated, however, as sug- 
gested by the reduced of the profile fit discussed in 
section 13.21 below. 

2.3. Photometric redshift measurements 

Of the 129 groups in our full sample, 95% con- 
tain two or more spectroscopically confirmed members, 
3% have one spectroscopically confirmed member, and 
the remainder have redshifts determined photometri- 
cally from the red seq uence of member galaxies (as in 
iFinoguenov eTalllMOil . The average redshift error for 
the group ensemble is ^ 0.0017, only slightly larger than 
their typical velocity dispersions of 300-450 km/s. 

Multicolor groun d-based imaging in over 30 bands 
(jCapak et al.l l2007l ) also provides photometric redshift 
information for all of the source galaxies along the same 
line of sight. We use photometric redshift estimates 
from the LePhare template- fitting code, which are 
updated from those published in lllbert et al.l ()2009[ ) by 
the addition of deep i/-band data, and small improve- 
ments in the template fitting technique. We have com- 
pared photo-zs from 10801 galaxies at z ^ 0.48, 696 
at z ~ 0.74, and 870 at z ~ 2.2 to spectroscopic red- 
shift measurements with the Very Large Teles cope (VLT) 
Visib le Multi-Object Spectrograph (VIMOS) (' Lillv et al.l 
l2007f ) and the Keck Deep Extragalactic Imaging Multi- 
Object Spectrograph (DEIMOS). The r.m.s. dispersion 
in the offset ctaz between photometric and spectroscopic 
redshift is 0.007(1 + z) at i\j^ < 22.5 and 0.02(1 -I- z) at 
i\g - 24 and z < 1.3 (or 0.06(1 + z) for i\g - 24 and 
z > 1.3). 

To mitigate against catastrophic failure in estimated 
photo-zs, for example due to confusion between the Ly- 
man and 4OOOA breaks, we reject from the sample all 
source galaxies with a secondary peak in the redshift 
probability distribution function (i.e. galaxi es where the 
param eter ZP_SEC is greater than zero in the lllbert et al.l 
120091 catalog). The rejected zp_SEC > galaxy popula- 
tion is expected to contain a large fraction of catastrophic 
error s (roug hly 40%-50% - lllbert et al.T2006l : lllbert etall 
l2009f ). For the purposes of cosmological constraints, we 
further exclude from the sample objects with relative red- 
shift uncertainties Az/(1 + z) > 0.05, taking the aver- 
age redshift error to be Az = (zu68_gal— zl68_gal)/2.0, 
where zu68_gal and zl68_gal are the 68% confidence 
limits on the redshift, ba.sed o n the photo-z probability 
distribution (jllbert et al.ll2009( ). Our final source sample 
consists of all galaxies passing these cuts that lie within 6' 
of a group center. Individual galaxies may enter into the 
final sample multiple times if they lie within 6' of more 
than one peak. The photo-z quality cuts reduce the num- 
ber density of source galaxies to 26 galaxies/arcminute^, 
for a total of 3.7 x 10^ galaxies (3.1 x 10^ in the re- 
stricted sample). The mean redshift of the final sam- 
ple is (z) = 0.95 and the mean relative error in red- 
shift is Az/(l-|-z) = 0.018, while the mean magnitude 
is (If814w) ^ 24. 
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3. METHOD 

3.1. The Stacking Method 

If we consider a source galaxy (or 'source' hereafter) at 
redshift zs being lensed by a foreground group (or 'lens' 
hereafter) at redshift zl and observed at zq, in the weak 
limit the tangential shear induced by the lens will be: 

jt{r) = [S(< r) - S(r)] /S, = AE(r)/S, (1) 

where !](< r) is the mean surface density interior to pro- 
jected (physical) radius r, S(r) is the azimuthally aver- 
aged surface density at r, and Ec is the critical surface 
density defined as: 



Ds 



AttG DlDls 



(2) 



Here Ds denotes the angular diameter distance from the 
observer to the source 



Ds = /fc(t^[2o, zs\)a{zs) , 



(3) 



where uj[zo, zs] is the comoving (or coordinate) distance 
along a radial ray between the observer and the source 



, r° cdt r 

Jzn Jas "'W Jas 



cda 



(4) 



Dl and Dls are angular diameter distances from the 
observer to the lens and from the lens to the source re- 
spectively, given by 

Dl^ fk{i^[zo,ZL])a{zL) , (5) 
DLS^fk{(^[zL,zs])a{zs), (6) 

(e.g. iBartelmann fc Schneiderll200lD . 

If we consider the case of a flat cosmology (fc = 0), 
/fc(a;) = uj and thus we can rewrite Ec more simply 
in terms of the comoving distances ujl = lu[zo,zl] and 
Us = uj[zo,zs]: 



UJS 



AttG uil{ojs ~ ul) 



[I + zl). 



(7) 



The critical density incorporates all the geometric depen- 
dence of lensing; it is the non-linear dependence of Ec on 
zs and zl that makes it hard to stack the signal from 
different source-lens pairs in a straightforward way. We 
can simplify the dependence, however, by defining the 
comoving distance ratio x = us/oj^. Note that x > 1 
for sources chosen to be behind the lens; we will consider 
objects with both x > \ and x <1 below. In terms of x, 



(1 + zl) 



1 



(8) 



47rG' UJL (1 - l/x) ■ 

the value of the 



We can also write this in terms of Ec,oo, 
critical density in the limit x oo : 



Ee = 



1 



where 



(l-l/x)^'°"' 

C2 (1 + Zi) 



47rG 



(9) 



(10) 



depends only on the lens properties, not on the properties 
of the source galaxy. 



From Eqn. [l] the geometry of all source-lens pairs now 
takes on a universal form 



7t(r)Ec^co/AE(r) = r(a;) = 1 - 



1 



(11) 



T{ x) corresponds, e.g ., to the lensing efficiency E defined 
bv lGolse et al.l (|2002| ). In as much as the measured tan- 
gential ellipticity St of each source galaxy is an estimator 
7t of the true tangential shear 74, we can construct a 
weighted sum of estimates from individual source galax- 
ies j with respect to lensing centers i to recover the uni- 
versal geometric dependence: 



u',j-7t,jj-S(c^oo),/AE,(r,j-) 



(12) 



with weights Wij chosen to maximize the signal-to-noise 
ratio or sensitivity to cosmological parameters, as dis- 
cussed below. 

Since we are just fitting the data to a fixed function, 
cosmology appears to have disappeared from Eqn. ((T^ . 
In fact, it is hidden in the conversion from measured 
source and lens redshifts to inferred source and lens dis- 
tances. For a given cosmology we convert redshifts to 
comoving or angular diameter distances, construct the 
weighted sum in Eqn. [121 and calculate the of r(a;) 
with respect to the theoretical expectation (1 — 1/a:). 
This gives us the relative likelihood of that particular set 
of cosmological parameters; iterating over this process 
then allows parameter constraints. The only remaining 
problems are to determine the optimal weights Wij and 
surface mass density contrast AE(r). We discuss these 
calculations in the next section. 

Before proceeding we should note that our simple 
stacking analysis ignores several complications. First, it 
ignor es the distincti on between true and reduced shear 
(e.g. IShapird l2009t ) . In the weak shear limit the two 
are identical, and for the groups considered here the sur- 
face mass density is low enough that the contribution 
from non-weak shear corrections is unimportant outside 
50-100 kpc. The effect of the second-order cor- 



rection term is illustrated in figure 4 of iLeauthaud et aTl 
()2010D . Its contribution is roughly comparable to that of 
the stellar mass in the central group galaxy (see Fig. [2] 
below) at large projected radii, and always less than the 
stellar contribution at radii less than ^ 50-100 h^^ kpc. 
Given only that ^2% of our sources lie at such small 
projected radii, the effect of these contributions on our 
fits should be negligible (excluding from the sample all 
sources within 15' of group centers^, for instance, has no 
effect on the final results). 

Second, we have also assumed fiatness in separating the 
dependence on the lens distance and the source distance. 
While current cosmological cons traints indicate an al- 
most completely flat universe (e.g. lLarson et al1l2011l find 
< 0.01 from various sets of constraints) , it would be 

^ Note that in order to evaluate Ax^ smoothly as we vary the 
cosmological parameters, the sample selection has to be indepen- 
dent of cosmology, and thus our cuts have to be in an observable 
such as angular separation, rather than a cosmology-dependent 
quantity such as physical separation. 
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nice to be able to relax this assumption. Unfortunately, 
while it is still possible to fit shear ratios between indi- 
vidual pairs in the general case, there is no simple way 
of stacking all measurements together into a single func- 
tional form, since the dependence on the two redshifts 
can no longer be factored out of fk{uj[zL, zs]) in a simple 
way. We can estimate the effect of curvature by consider- 
ing the series expansion for fk in the limit <C 1. For 
typical values of cj ~ c/iJg, the next term in the series 
is smaller by a factor l/6(a;/i?o)^ ^ 1-6 x 10""^, where 
Rq = c/(iJovK^)- Thus in realistic non-flat cosmolo- 
gies, we expect a correction of order 1.6 x 10^^ to our 
values of r(a;). Compared to our cosmological sensitivity 
Ar ~ 0.05 (see section |4?2] below) . this represents a 3% 
correction to our derived parameters. This correction is 
smaller than errors on fix we obtain below, although it 
would be quickly become important in larger surveys. In 
what follows we will ignore the complication of non-zero 
spatial curvature. 

Third and lastly, our analysis assumes a specific func- 
tional form for the surface mass density contrast AI](r), 
namely a projected NFW profile. We will show below 
that this functional form is in fact an excellent fit to the 
stacked data. We could use instead an empirical profile 
based directly on the data itself, but given the agreement 
between the NFW model and the data, this would not 
affect our results significantly. 

3.2. Radial Dependence and Optimal Weighting 

Our goal is to measure the redshift dependence of 
the group lensing signal. The radial dependence of the 
surface mass density contrast AE(r) around groups, al- 
though intrinsically very interesting, is essentially a nui- 
sance parameter in this calculation. We need to deter- 
mine AI](r), however, in order to weight measurements 
of individual source galaxy shapes optimally when esti- 
mating r(a;). The density contrast a r ound groups was 
studied in detail in iLeauthaud et al.l (|2010f ) ; we repro- 
duce the same calculation here, stacking with respect 
to physical radius the signal from all groups with well- 
determined redshifts and centers. 

The density contrast profile includes contributions 
from four main terms^: the weak shear contribution of 
the main halo (the '1-halo' term), the average weak shear 
contribution from nearby halos (the '2-halo' term), a 
weak shear contribution from the stellar mass of the cen- 
tral group galaxy, and the second-order corrections to 
the shear in the center of the main halo. Of these, only 
the 1-halo term is important here; the 2-halo term only 
becomes significant at large projected radii (r ~ 4 hj^ 
Mpc), while the stellar and second-order terms are only 
significant at small projected radii (r < 50-100 h^^ kpc) 
where we have very few galaxies in the source sample, as 
discussed in Section [Ol 

We expect that the 1-halo term for a single group 
should follow a projected Navarro-Frenk- White (NFW) 
profile S(r), whose form f n fw (t I ''' s) is given, e.g. in 
IWright fc BrainerdI (|2000l ). The profile has two free pa- 

^ In some cases, a fifth component might be necessary to ac- 
count for the mis-identication of the central galaxy l|Johnston et al.l 



rameters, a scale radius and a normalization Sq, or al- 
ternately it can be defined in terms of a virial radius r^r 
and a concentration c = Vmrl'^s- The expected values of 
these parameters can be inferred from X-ray fiuxes, X-ray 
scaling relations, and theoretical concentrati on-mass re- 
lation s. Using the co ncentration relati o ns of jZhao et al.l 
(|2009f ). for instance, ILeauthaud et al.l ()2010D predicted 
concentrations in the range 3.6-4.6 for the COSMOS 



R 



groups (these values assume the definition 

thus they correspond to the values C2ooc from'Zha o et al.l 
|2009( ). Simulations suggest that individual halos will 
have significant (^50% ) scatter around these mean val- 
ues (jZhao et al.ll2009f ). Finally, we note that here we 
are considering an average profile for all groups in the 
sample, where the averaging is weighted by surface mass 
density contrast AE. At a fixed redshift and fixed con- 
centration, Ts and r-uir for each group should scale as 
Ml/3 and our stacked profile would be similar to a mass- 
weighted average. Variations in concentration and red- 
shift complicate this behavior, but we can still use mass- 
weighted averaging to guide our expectations as to the 
final values for the concentra tion or scale radiu s. Using 
the concentration relations of IZhao et al.l ()20Q9l) and as- 
suming a WMAP7 cosmology, for instance, we predicted 
a mass- weighted average scale radius of = 154 K^^ kpc 
for our groups. 

The surface mass density contrast AI](r) is related to 
the tangential shear by 



AI](r) =E, X7t(r). 



(13) 



Thus AS](r) can be estimated as in Eqn. 8 of 
ILeauthaud et al.l (|2010f) 



AE(r) 



(14) 



[20071) . In our case, we neglect this term which is sub-dominant in 
our sample because we have optimized the centering using weak 
lensing (George et al in prep.). 



Fig. |2| shows the surface mass density contrast for our 
stacked sample of groups. Points with error bars in- 
dicate the mean value inferred from Eqn. 1141 binned 
logarithmically in radius. The thin dotted lines indi- 
cate contributions from a projected NFW profile (middle 
blue line), the stellar mass of the central galaxy in the 
group (lower red line), and the sum of these two com- 
ponents (upper green line). The normalization of the 
stellar contribution is based on the mean stellar mass of 
the central group galaxy, as inferred from photometry 
(see ILeauthaud et al.ll2011l for details). 

The main part of the profile is well fit by a pro- 
jected Navarro-Frenk- White (NFW) profile E(r) with a 
scale radius ~ 160-180 h^^ kpc and a normalization 
So ~ 200 h^Q Mq pc~^ at rj. The best fit values of these 
two parameters are strongly correlated, so we choose in- 
stead to normalize the profile at a radius rfix — 250 h^^ 
kpc where the measurement errors are small and the am- 
plitude is less dependent on r^. Thus the profile can be 
written 

^[r)==^o- -, j-^. (15) 

This gives us constraints on Sq — S(rfix) and which 
are more independent of each other. 

We determine values for Sq and Vg by calculating 
with respect to a projected NFW model. We do not bin 
the data in radius, since the bin boundaries would change 
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Fig. 2. — The surface mass density contrast as a function of 
projected physical separation, stacking all groups with well-defined 
redshifts and centers (the 'restricted' sample). The thin dotted 
lines indicate contributions from a projected NFW profile (middle 
blue curve), the average stellar mass of the central galaxy in the 
group (lower red line), and the sum of these two components (upper 
green curve). 

with cosmology or scaling, but sum the contribution from 
each galaxy individuahy. For the full group sample and 
assuming WMAP7 cosmology, we obtain the lowest re- 
duced by considering all sources at x > 1.1, with no 
other cut on redshift errors. This gives a slightly high 
value for the scale radius — 187^29 ^70 ^P^^ however. 
Considering only the restricted sample and/or sources 
with Az/(1 + z) < 0.05, we obtain values of closer 
to the expected value. This suggests that the remaining 
20% of the group sample may be affected by centering 
problems, or that it may include many unrelaxed or dis- 
turbed groups with lower mean concentrations. The de- 
pendence on photo-z cuts could indicate that the lensing 
signal in all groups is diluted by group member contami- 
nation in the source population when the limits on photo- 
z errors are relaxed. For the restricted group sample with 
the cut on source redshift errors, we obtain best fit values 
Eo = 98.8 ± II/17OM0PC-2 and = 158tg h^^ kpc. 
The reduced is marginally higher for this sample than 
for the much larger uncut sample, but the best-fit value 
of Ts closer to the expected value, so we will take this 
as our fiducial profile, and marginalize over values of 
and So in this range for our cosmological calculations. 
The best fit value of also places some constraints on 
possible centering errors for the groups. We have tested 
the effect of centering errors by adding random offsets 
to the individual group centers, with r.m.s. values of 6", 
12", and 24" in each coordinate. The resulting profiles 
are still well fit by our model, but the best-fit value of 
Ts increases to 220, 260 and 340 h^Q kpc for the three 
cases respectively. This suggests average centering errors 
are < 6" ~ 25-50 h^Q kpc in each coordinate, consistent 
with other estimates of the centering uncertainty (George 
et al. in prep.). 

We have also investigated other forms of stacking. In 
principle we could correct for the predicted variations 



in concentration, for instance, stacking in r/r^, or we 
could stack in comoving rather than physical coordinates. 
Testing scaling the radius by r^, r^r, or (1 -I- z^), we find 
little or no significant improvement in the of th^ fit to 
the radial profile. In particular, we find only a marginal 
indication of any trend in concentration with mass or 
redshift. Given that the halo-to-halo scatter measured in 
simulations is comparable or larger than to the average 
trends over the mass and redshift range spanned by our 
group sample, this is perhaps unsurprising. Furthermore, 
since the concentration relations are themselves depen- 
dent on cosmology, we would have to account for this 
dependence in our marginalization over cosmological pa- 
rameters, so we will not attempt to correct for variations 
in concentration. We can also consider other analytic 
fits to the profile. We find that NFW is preferred over a 
cored isothermal profile at the 95% confidence level, and 
preferred over a singular isothermal (E(r) cx r~^) profile 
at 97-98% confidence. Thus, our stacked profile provides 
empirical confirmation of the NFW model, in agreement 
with other recent h igh-precision measu rements of clus- 
ter density profiles (jUmetsu et al.|[2blll Okabe et al. in 
prep.). 

Finally, we note that the fit to the radial profile gives 
us an independent check of our empirical shear variance 
estimates. Because we fit the profile without binning, we 
have very large number of degrees of freedom and thus 
a narrow range of expected scatter in the reduced x^- 
The best-fit NFW profile has a reduced x^ of 0-931 with 
a very small (0.0027) expected scatter, so we conclude 
that our empirical variance is probably over-estimated 
by ^ 7%, corresponding to error bars which are 3.6% too 
big. We correct for this in all our subsequent analysis, 
multiplying the empirical shear variance by a factor of 
0.931. 

4. RESULTS 

4.1. The Geometric Signal 

Given a functional form for the radial dependence of 
the surface mass density contrast AE(r), we can proceed 
to estimate r(x) via Eqn. [121 The weights in the sum 
can be calculated as the inverse variance of the Fy 

w,, = (var[r,,])-i ^ ("^plMV var[7M,r' • (16) 

where the variance of the tangential shear is determined 
empirically, as described in section [2.21 

The 'model' here, the geometric sensitivity function 
(1 — 1/a;), is fixed, while the data vary as we change Sq 
and r^, which both change AI](r), and the cosmolog- 
ical parameters, which map the redshifts (2:5,2;^) onto 
a;- values and also determine Sc,cx)- If we restrict our- 
selves to flat cosmologies with two components, matter 
and dark energy, then the goodness-of-fit depends on Sq, 
Ts, r^x, and w. 

Fig. [3] shows r(a:;) over the range x — [0, 5] for a 
WMAP7 cosmology with = 0.27, I^a = 0.73. The 
points are weighted averages in bins of 0.3 in x, while the 
solid (red) curve is the theoretical expectation: 

r(x)=0 for x<\\ (17) 

= 1 — 1/a; for a; > 1 . 

Weights here are inverse variance, as in Eqn. 1161 The 
errorbars on the data points are calculated as usual for 
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Fig. 3. — The geometric dependence of the lensing signal, plotted 
versus the distance ratio x = uj^/^L^ for ^ WMAP7 cosmology 
with f2jvf = 0.27, S^A = 0.73. The solid (red) curve shows the 
theoretical prediction r(x) = for a; < 1; r(x) = 1 — 1/x for 
a; > 1. 

an inverse- variance weighted average: 

OT.bin = ^ Wi] , (18) 

Vbin / 

where the sum is over ah pairs with values of x in 
the bin. 

Clearly the geometric signal is present in the COS- 
MOS data, and measured to reasonable significance over 
a broad range of distance ratio x. Given the possible 
systematics in the measurement discussed below, the ex- 
cellent agreement between theory and data illustrates the 
potential of the method. On the other hand, "x^-by-eye" 
is somewhat misleading for this figure, as the binning 
in X may hide systematics at particular distance ratios. 
(There is a 1.9a indication of positive signal in one bin 
sd X < 1, for instance. This could indicate photo-z errors 
are scattering sources to lower redshifts, but given the 
number of bins there is a 40% chance that this is simply 
a random statistical fluctuation.) As discussed below, 
in a small field structures at a few redshifts can domi- 
nate the lensing signal, introducing excesses or deficits of 
mass along the line of sight that dominate the signal at 
particular values of x. The bins chosen here are broad 
enough to smooth out many of these features, but clearly 
a goodness-of-fit measurement over the whole data set is 
required to determine the statistical significance of the 
apparent agreement in Fig. 3. 

4.2. Parameter Constraints 
We can estimate goodness-of-fit by calculating 

X^=^mj[{l~l/x,,)~T,,]^ (19) 

This sum depends on Sqj and the cosmological pa- 
rameters, so marginalizing over the first two parameters 




12 3 



Fig. 4. — The sensitivity function Ar{zg , Z]^). Points with high 
source or lens redshifts are the most sensitive to changes in cos- 
mology. 

gives constraints on the dark energy density iljc and the 
equation-of-state parameter w. While the resulting value 
of will tell us whether the data are a good fit to the 
model, this is not necessarily the most sensitive way of 
determining cosmological parameters. In particular, for 
fiat (fc = 0) cosmological models with a cosmological con- 
stant with density parameter fix — ^A, as the value of 
JIa increases all distances will increase, and therefore so 
will X and (1 — 1 /x) . Thus the signal from a given source- 
lens redshift pair will be compared to r(x) at a value of x 
which is larger by some factor. The most sensitive probes 
of this re-scaling will be points at large x. To constrain 
ilx more precisely, individual measurements should be 
weighted by this sensitivity. 

The exact sensitivity to cosmology itself depends on 
the cosmological parameters. We can estimate a sensi- 
tivity factor and thus a weighting that will be close to 
optimal over the whole range of JIa, however, by calcu- 
lating 

Ar(zs, zl) = Tqi [x{zs,zl)] - Tog [x{zs, zl)] 

= [xo9{zs,zl)V^ - [xqi{zs,zl)]^'^ (20) 

where Fqi Fpg are the model F(a;) evaluated for cosmolo- 
gies with f^A = 0.1 and 0.9 respectively (with w = —1 in 
each case). Fig. |4] shows this sensitivity function ver- 
sus (zs,2l)- (Note that this sensitivity f unction has 
been calculated previously, e.g. in Figure 3 of lGolse et al.l 
120021 ). For zs < z^, the sensitivity function is zero since 
the model value T{x) = independent of cosmology. For 
zg > Zl, sensitivity generally increases with source or 
lens redshift. 

To maximize signal-to-noise in our cosmological con- 
straints, we apply this weighting quadratically in our 
previous expression for x^' 

ij ij 

This can be converted to a likelihood by assuming 
the error distribution is Gaussian. This is only approx- 
imately true in our case, but determining more accu- 
rate error distributions would require a significantly more 
complex error analysis, so we will leave this to future 
work. Fig. [5] shows the likelihood function for JIa, calcu- 
lated using the restricted group sample and normalized 
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Fig. 5. — Likelihood function for with a uniform prior on Ts 
over the range 120-200 h-^ kpc, using the restricted group sample. 
(The likelihood has been normalized so the area under the curve is 
1.) The dashed vertical line indicates the value Q/^ = 0.848 where 
the likelihood peaks. Dotted vertical lines indicate 68.2, 95.4 and 
99.7% (1, 2, 3-cr) confidence regions. The solid vertical line and 
shading indicate the mean W MAP seven-year valu e = 0.727 
and 68% confidence range from lLarson et aLl 120111 ). The vertical 
green line indicates the value fi^ = 1/3, the lower limit for cosmic 
acceleration (qo < 0). 

SO the area under the curve is 1. The equation of state is 
fixed to w — —1 and we have marginahzed over Sq and 
Vs. The dashed vertical hue indicates the value where 
the likelihood peaks. Dotted vertical lines indicate 68.2, 
95.4 and 99.7% (1, 2, and 3-a) confidence regions. The 
solid vertical line and shading indicate the mean WMAP 
seve n-year value r ^ A = -727 and 68% confidence range 
from lLarson et all ()20lH ). 

The results are insensitive to the priors in Sq provided 
they are reasonably broad, but they depend strongly on 
the smallest values of considered. This is because 
of the dependence of inferred physical distance on cos- 
mology. For low values of fix, the transverse distance 
inferred from a given angular separation on the sky is 
smaller. Conversely, if we fix Ts to a small value, small 
values of fix are preferred. Because the radial varia- 
tion of the profile over the range of our data is much 
stronger than the redshift variation, the radial fit drives 
the values and thus biases our cosmological results if 
unphysically small values of are allowed. The aver- 
age scale radius is predicted to be ~ 160 h^Q kpc for 
our sample, based on theoretical concentration relations. 
Simulations sho w ^50% scatter i n concentration from 
halo to halo (e.g. IZhao et al.lf2009t iReed et al.l[201lL and 
references therein), but most of these variations should 
average out in the set of ^100 objects considered here, 
provided they represent a reasonably unbiased sample. 
On the other hand the conversion from angular to phys- 
ical radius will vary by ~10-15% for the range of cos- 
mologies considered, and the range of uncertainty on our 
fitted value of is ~ 130-210 /ifo^ kpc. Thus, we take a 
uniform prior on over the range 120-200 h^^ kpc, al- 
lowing for a variation of ±25% around the fiducial value 



Ts — 160 h^Q kpc. Extending the range of our priors to 
higher values of has little effect on the results, while 
extending it to lower values of decreases the lower limit 
on fix- 

Overall, we obtain the estimate flx{w — —1) = 
J7a = 0.848^q'587^, corresponding to a detection of 
dark energy at more than 99% confidence. This value 
is consistent with the most recent WMAP analysis of 
CMB anisotropics, which finds a mean value 0.727t|^:j;^!^ 
(jLarson et all 1201 iD . Our 68% confidence range is ap- 
proximately 4 times wider than that of WMAP; given 
the small size of the field considered here (1.64 deg^), 
however, this level of precision demonstrates the power 
of the geometric test. We note however that our error 
estimates do not include systematic effects. We estimate 
the magnitude of some of the possible systematics in the 
next section. 

In a cosmology with multiple components with 
equation-of-state parameters w and density parameters 
fitu , cosmic dynamics can be characterized by the decel- 
eration parameter 



(22) 



For a flat universe with components flm and JIa, we find 
a value qo ~ ^0.77lo.o66 the present day. Present-day 
acceleration, which corresponds to go < or fl^ > 1/3 if 
w = —1, is detected at greater than the 98% C.L. (solid 
[green] vertical line in Fig. O. 

We can also extend the constraints to more general 
dark energy models with w —1. Fig. |6] shows the like- 
lihood surface for models with a dark energy component 
fix with an equation-of-state parameter w. Our cur- 
rent results provide an upper bound oi w < —0.4 (68.2% 
C.L.). They do not provide a lower bound, although 
there is some information in the constraints on w as a 
function of flA- The shape of the confidence regions is 
similar to t hose derived by g eometric tests using strong 
lensing fe.g. lJullo et al.l [2"010l), although our contours are 
shifted upwards to less negative values of w, perhaps be- 
cause of the redshift distribution of our lenses. The shape 
of the confidence regions also di ffers from those de rived 
from obser vations of the CMB (jLarson "erall[201ll). su- 
perno vae (Sullivan et al, 1201 If ), or BAO (Per cival et al.l 
l2010t [Beutler ct_al, 201 ll ) , providing interesting comple- 
mentarity with these other methods. 

4.3. Cosmic Sample Variance and Systematics 

4.3.1. Basic Noise Sources 

There are several basic sources of noise in our mea- 
surement of r(a;). Individual galaxy shapes are noisy 
estimators of the tangential shear, and shape measure- 
ments themselves are not perfect, particularly for small 
or faint galaxies. The redshift errors in our lens popula- 
tion are negligible (~ 0.0017 on average), but the pho- 
tometric redshift errors for the source population may 
contribute significantly to the errors, once again particu- 
larly for faint galaxies at high redshift. These statistical 
errors are included, at least approximately, in our esti- 
mates of x^- likelihood estimator also assumes a 
particular form for the average radial density profile of 
the groups Eq. This profile may not correspond exactly 
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Fig. 6. — Likelihood surface for the cosmological parameters Qx 
and w (grayscale). The dotted contours indicate regions enclosing 
68.2%, 95.4% and 99.7% (1, 2, S-cr) of the probability in the region 
covered by our priors (—2.5 < ui < —0.1, < fix < 1). The "X" 
marks the mean WMAP seven-year value from lLarson et aP l)2011h 

to the true mean surface mass density contrast, and indi- 
vidual groups may deviate from the average in ways that 
depend systematically on redshift. The tests described 
in section 13.21 suggest these effects are at or below the 
level of our statistical errors. 

The relative importance of the source or lens sample 
selection is less clear a priori. To explore the depen- 
dence of our results on source galaxy properties we have 
rerun parameter constraints with samples cut in magni- 
tude range, S/N, size, and redshift error. In each case 
reducing the sample size increases the final error in n^, 
but there is no obvious systematic trend for cuts in mag- 
nitude, S/N or size, beyond the dependence on source 
redshift discussed in section 14.3.21 below. The precision 
of the final results is quite dependent on the redshift 
cuts we make in the source sample and on cuts to the 
group sample; including all the groups in the full (rather 
than the restricted) sample increases the uncertainties 
by ^-^25%, including source galaxies with zp_SEC> 
increases them by '^30%, and including sources with 
Az/(1 + z) > 0.05 doubles them. Examining the cat- 
alogue, it seems the particular sensitivity to redshift er- 
rors is related to sources at z — 2-3, which contribute 
strongly to the constraints due to the sensitivity func- 
tion (cf. Fig. in . Cutting these objects out, or diluting 
their signal with many poorly determined redshifts, both 
reduce out final sensitivity quite strongly. We have also 
experimented with using more restrictive cuts on redshift 
errors, or using a larger aperture 10' aperture around 
each group; these do not improve our results significantly. 

Finally, the relative importance of redshift errors is also 
slightly unclear in our method. We have attempted to 
estimate the contribution of photometric redshift errors 
to the total error budget by adding scatter to the photo- 
metric redshifts in the source catalogue. Adding scatter 
equal to 5x or lOx Az (increasing (Az/(l-|-z)) to 0.09 or 
0.18) increases our 68% C.L. bounds on 51a by 12% and 
30% respectively. This suggests that our current redshift 
errors contribute only 5-10% of the total error budget. 



If so, this is a promising result for other surveys which 
are likely to have redshift errors larger than COSMOS. 

4.3.2. Cosmic Sample Variance - Empirical Estimates 

A more complicated source of systematic uncertainty 
in the measurement comes from structure along the line 
of sight. The COSMOS survey probes a relatively small 
field, and several large structures are clearly visible in 
the galaxy redshift distribution below redshift z = 1. We 
can look for evidence of cosmic sample variance effects in 
various ways. One simple estimate of the importance of 
cosmic variance can be obtained by boot-strapping, split- 
ting the sample into two or more disjoint sets. We can 
test for systematics by measuring how Agg , the width of 
the 68% confidence region for Ha , changes relative to the 
Poissonian expectation Ags c>c i/nT/rio when we cut the 
sample from tiq sources down to ric- Splitting the field 
into 4 quadrants with roughly equal numbers of source 
galaxies, we find Agg increases by a factor of 1.9-2.5, so 
this seems consistent with the factor -y/ric/no = 2 ex- 
pected from Poisson scaling. We note however that one 
quadrant gives a best fit value of J7a = (albeit with a 
68% uncertainty of -1-0.6) whereas the others give values 
of ~0.9. Examining the lenses and source distribution 
in this quadrant, it is not immediately obvious whether 
specific structures produce this shift. This is a sobering 
lesson about the possible effects of cosmic sample vari- 
ance. 

Splitting our lenses into two groups by redshift, each 
with Uc ~ ?^o/2, we find Agg increases by a factor of 1.33 
for zl > 0.4 and 2.4 for < 0.4. Thus the noise in- 
creases at a roughly Poissonian rate in the high redshift 
sample, while in the low redshift sample it increases much 
more quickly. Similarly in a sample cut at zs > 1, Agg 
increases by 1.33, while for 2:5 < 1, it increases by 1.73. 
These results suggest the signal at low redshift is more 
prone to systematics. We have tested alternative weight- 
ing schemes that attempt to correct for the trend with 
source redshift. Down-weighting sources with zs < 1-5 
by a factor of 0.5 or 0.33 reduces Agg by 30% or 40% re- 
spectively, but moves the peak of the likelihood 1.3cr or 
1.7a away from the WMAP 7 value. Thus, it seems there 
is some trade-off between precision and accuracy in the 
redshift weighting. Clearly this subject requires further 
theoretical work, using realistic simulations of large-scale 
structure. For the moment, in the absence of an optimal 
weighting scheme motivated by theory, we choose not to 
apply either weighting to our final results. 

One final concern is that the average line of sight to all 
our groups could be slightly over- or underdense. This 
would introduce a baseline shift in our model of r(j:), 
modifying it to: 

T{x) = Fo for a; < 1 ; 

= (1 - 1/x) + Tq for x>l. (23) 

We can test for a constant offset Fg ^ most easily in 
the range x = [0, 1] where the expected signal is zero. 
We find Fq = 0.02 ± 0.03, so there is no significant evi- 
dence for an offset. Furthermore, adding Fq ~ ±0.02 to 
our model gives essentially identical constraints on D,\ 
(the peak value shifts by less than O.lcr), so the method 
appears to be robust to any small offset of this kind. 
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Fig. 7. — The predicted contribution to cosmic shear errors 
from unidentified structure along the line of sight up to a given 
source redshift. The contribution is integrated over a 6' aperture, 
assuming a WMAP7 cosmology and the COSMOS source redshift 
distribution. 

4.3.3. Cosmic Sample Variance - Theoretical Prediction 

We can also use the error description of iTavlor et al.l 
(|2007f ) to estimate the effect on our measurements of 
structure along the line of sight. They calculate that the 
tangential shear induced by large-scale structure between 
the observer and two background galaxies at redshifts Zi 
and Zj introduces a covariance in shear measurements 
given by 

(24) 

where we have integrated over a circular aperture of ra- 
dius 9 by multiplying by the Fourier transform of the 
aperture (the term in brackets), and C ]^ (£ ) is the tomo- 
graphic cosmic shear power spectrum ()Hulll999l) . 

In our case, we are only concerned with the auto- 
correlations between redshift bins with i = j. These 
will give an estimate of the excess variance added to our 
shear measurements by cosmic structure, as a function 
of source redshift. Fig. [7] shows this error term calcu- 
lated for an aperture of 6' in a WMAP7 cosmology and 
assuming the COSMOS redshift distribution. We have 
calculated the shear error using 20 discrete bins in red- 
shift between < z < 2. Ideally we would use a contin- 
uous cosmic shear in this measurement, as described in 
iKitching etahl (|201lD . but since the error contribution 
we find here is small and smoothly varying, this approx- 
imation seems adequate. 

We see that the extra error term is always less than 
our empirical shear dispersion > 0.25, and that it 
reaches a maximum of '~6% of the empirical dispersion. 
This suggests the contribution from cosmic shear is much 
smaller than the excess variance seen in the previous 
section, which may then be due to individual halos or 
to other systematics. Clearly more detailed simulations 
are needed to determine realistic cosmic sample variance 
errors for our particular technique, but these estimates 
reassure us that systematics do not completely dominate 
our current results. 

5. CONCLUSIONS 

Gravitational lensing is a versatile source of cosmologi- 
cal tests; it can be used to measure both the spectrum of 
matter fluctuations in the universe, and also its geometry. 



One particularly simple geometric test measures the rate 
at which the signal grows behind a lens and uses this to 
determine empirically the relationship between angular 
diameter distance and redshift. Here we introduce a new 
form of this test, the 'stacked shear ratio test', which can 
be applied to large sets of individually weak lenses. By 
factoring out the source and redshift dependence from 
the expected amplitude of the lensing signal, we define a 
universal geometric dependence T{x) = (1 — 1/x), where 
X = ujs/ujl and ujs,ujl are the comoving distances to 
the source and lens respectively, which can be estimated 
from individual measurements of tangential shear. As 
we vary cosmological parameters, the goodness-of-fit of 
our data with respect to the universal form provides a 
constraint on the geometry or equation of state of the 
universe. Physically, the test measures the relative dis- 
tance scale behind the lenses; a higher value of fix will 
increase this scale, so we can get particularly sensitive 
upper bounds on the density of dark energy. 

The new method requires a significant lensing signal, 
although the signal can be spread over many individ- 
ual lenses, and a large set of sources with well deter- 
mined redshifts. The X-ray-selected group sample from 
the COSMOS survey provides a good test of the method, 
although cosmic sample variance is a concern over a field 
this small. Stacking groups with well-determined cen- 
ters, we find a radial surface mass density contrast in 
good agreement with theoretical predictions. Examin- 
ing the mean tangential shear behind the stacked group 
sample as a function of x, we see a clear manifestation 
of the geometric signal. 

Calculating goodness-of-fit for flat cosmological mod- 
els with matter and a cosmological constant, with opti- 
mal weighting for sensitivity to cosmology over a broad 
range of fix, we obtain the constraints fl\ = flxiw = 
-1) = 0.8481^ ?^^ (68.2% CL). Extended to dark en- 
ergy models with constant equation-of-state parameters 
w ^ —1, our results do not constrain w from below, but 
give the upper bound w < —0.4 (68.2% CL). The shape 
of the confidence regions differs from those derived from 
observations of the CMB, supernovae or BAO, providing 
interesting complementarity with these other methods. 
Systematic uncertainties and the effect of cosmic sample 
variance are hard to estimate accurately without more 
detailed simulations, but the simple tests presented in 
Section 14.31 suggest they do not dominate our statisti- 
cal errors. Overall, the consistency of our results with 
more accurate determinations of VIa ([Larson et al.ll201l[ ) 
provides a good demonstration of the potential of the 
stacked shear ratio test in larger-area surveys. 

Our method does have several limitations. The sim- 
ple scaling only applies for flat cosmologies; for fc ^ 
the geometric term in the lensing equation depends sep- 
arately on zs and zl, so a more complicated 2-D fit is 
required. We have also worked exclusively in the weak 
shear limit where k ~ 7 4; 1. This does not restrict us 
significantly for COSMOS groups, since they are in the 
weak regime well into their central regions where other 
effects (e.g. baryonic mass, centering errors) become im- 
portant anyway. Applying our method to more massive 
clusters would require modeling the full shear, and thus 
introduce separate dependence on zs and zl- We also 
assume a specific functional form for the surface mass 
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density contrast, namely the projected NFW profile. In 
practice this form seems to be an excellent fit to the 
data, however, and therefore we would not expect our 
results to change significantly if we used instead of this 
an empirical profile based on the data alone. Finally, our 
lower bound on Qx is sensitive to the average physical 
scale radius of the lenses; sensible priors and/or reduced 
weighting at small radii are necessary here to avoid un- 
physical solutions which degrade the lower bound on fix ■ 
For the large samples of relatively low-mass lenses 
such as those expected from current (e.g. CFHTLenS^) 
and forthcoming lensing sur veys (e.g. DES^; Pa nStarrs^; 
LSST^; HSC^ Euclid^ see iMassev erall[20inl for a fuU 
listing), our method provides a simple, intuitive measure 
of geometry. The statistical uncertainties on J^a obtained 
for the 1.64 deg^ COSMOS field are only 4 x larger than 
those obtained from current CMB results ()Larson et al.l 
1201 It ), and should scale as (survey area)^/^. Thus if 
larger shear surveys can overcome systematics, there are 
good prospects for significantly reducing the uncertain- 
ties on abundance of dark energy. Given its sensitivity 
to sources at high redshift, our method might profitably 
be a dapted to lensing measurements bas ed on magnifica- 
tion (IVan W aorbokc ct al. 2010; Hildcbr andt et al.ll201lt 
iHeavens fc^ Joachimil 120110 . Since these require only 
measurements of magnitudes, not shapes, they may pro- 
vide larger samples of galaxies at z = 2-3. It could 
also be combined w ith CMB lensing (iHu et al.l 120071: 
iDas fc Spergell I2OO90 . which would provide an anchor 
point at the largest observable values oi x. In the longer 
term, as larger data sets become available and the sen- 
sitivity of lensing tests increases, their most interesting 
cosmological application might be to search for changes 
in the equation of state at particular redshifts. From this 
point of view, the method proposed here is flexible, as the 
redshift sensitivity depends simply on the distribution of 
source and lens redshifts used, and thus it can provide 
constraints over a range of redshifts, depending on the 
samples used. 
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